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A  computer  program  was  written  that  simplifies  the  design  of  parametric 
sonars.  This  program  has  sinceX^een  modified  tio  update  the  calculation  for 
secondary  directivity  index  and  to  include  calculations  for  secondary  beam- 
width  and  primary  beamwidth.  The  tolerance  onjthe  determination  of  secondary 
source  level  has  been  reduced  fromjjl.O  dB  to03.5  dB  with  a  consequent 
tightening  of  tolerances  on  all  other  related  calculated  quantities.  Finally, 
a  new  plotting  method  has  been  introduced  to  give  plots  in  the  form  of 
families  of  curves  for  ease  of  design  analysis. 
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INTRODUCTION 

A  computer  program  aimed  at  automating  the  design  process  for  para¬ 
metric  sonars  was  writtenJ  The  program  allows  the  designer  to  work  from 
a  known  secondary  source  level  (LSS)  and  a  known  secondary  frequency  (F) 
for  a  range  of  values  of  projector  area  (A),  primary  source  level  per  tone 
(LSP)  and  a  given  stepdown  ratio  ( FO/F) .  The  program  will  build  a  matrix 
of  possible  designs  that  can  be  compared  to  other  factors  to  achieve  a 
workable  and  realistic  design.  This  computer  program  is  based  on  the 
Mellen  and  Moffett  theory. 2 

Some  new  information  in  the  area  of  directivity  index,  a  tightening 
of  the  tolerances  on  the  computed  data,  and  an  alternate  routine  for 
plotting  data  have  been  incorporated  into  the  original  program.  Since  the 
computer  program  and  its  use  has  been  discussed  in  detail  in  the  earlier 
memorandum,'  only  the  changes  and  their  effects  will  be  taken  up  in  this 
memorandum. 


CHANGES 

The  program  remains  basically  the  same  as  that  discussed  in  the 
earlier  memorandum.'  The  changes  are  a  result  of  new  information  about 
the  directivity,  a  desire  to  tighten  the  tolerances  on  the  computed  informa¬ 
tion,  and  a  need  for  a  plot  that  arranges  the  data  in  the  form  of  a  family 
of  curves  for  ease  of  analysis. 

Directivity  Index  and  3  dB  bandwidth 

Investigation  has  revealed  that  the  original  expression  for  directivity 
index  by  Mellen  and  Moffett2  does  not  always  predict  the  secondary  frequency 
directivity  index  (NDIS)  accurately.  A  new  expression  was  developed  by 
Moffett^  which  was  deemed  more  accurate  over  the  region  of  concern.  This 
expression  is 

2 

NDIS  =  NDIP  -  10  LOG  (ye”°-65V(F/F0)  +1/2)  (i) 

where  Y  *  ttxF/C 5(GL) ( F0) ]  +  0.685  (AL)(R0)(F0)/F  (2) 

F  =  Secondary  Frequency  in  kHz^ 

F0  =  Primary  Frequency  in  kHz^ 

NDIP  =  Primary  directivity  index  in  dB1 
GL  =  Parametric  gain  in  linear  units 
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x  =  Primary  saturation  number.  In  this  program  a 
table  of  these  numbers  is  created  for  use.1 

AL  =  Absorption  loss  in  nepers/meter ^ 

A  related  term  that  is  now  calculable  is  the  secondary  frequency  3  dB 
beamwidth  (BWS).  This  is  expressed  as 

BWS  =  BWP  ^Ye'0*®5Y(F0/F)Z+l/2  deg.  (3) 

As  can  be  seen  from  the  expression,  BWS  depends  on  the  primary  frequency 
3  dB  beamwidth  (BWP).  In  this  program,  BWP  is  selected  as  being  for  a 
circular  piston  and  this  is  expressed  as _ 

BWP  =  (3.24)  (5.73)  jJ 4ttA( FO ) 2 ( 1 03 ) 2/C^  deg.  (4) 

where  C  =  Sound  velocity  in  meters/second. 

The  expressions  for  NDIS  and  BWS  contain  an  aperture  correction  orig¬ 
inally  derived  from  circular  piston  theory.  -However,  the  expressions  are 
also  applicable  to  square  piston  projectors.  Thus  all  of  the  computed 
information  relating  to  beamwidths  and  NDIS  applies  to  circular  and  square 
piston  projectors. 

A  complete  computer  printout  of  a  sample  example  is  shown  in  Appendix 
A.  The  sample  example  is  the  same  one  used  in  reference  1  so  that  the 
difference  between  the  new  and  old  programs  can  be  compared.  The  NDIS 
column  in  the  output  tables  now  reflects  the  new  theory.  Also,  two  new 
columns  appear,  one  for  BWS  and  the  other  for  BWP. 

Tolerance 


In  order  to  solve  a  given  problem,  the  computer  program  must  compare 
each  calculated  value  of  LSS  for  each  data  point  against  a  nominal  secondary 
frequency  source  level  (LSSI)  value  that  is  put  into  the  computer  in  the 
input  data.  The  tolerance  on  this  comparison  was  +1.0  dB  in  the  old  form 
of  the  program.1  The  new  tolerance  is  +0.5  dB.  TFus,  the  tolerance  on  all 
other  computed  values  that  depend  on  LSF  is  proportionately  tighter  than 
under  previous  conditions.  These  quantities  and  their  new  tolerances  are: 

G  =  Parametric  Gain  ^0.5  dB, 

LSP  =  Primary  Frequency  Source  Level  per  tone  +0.5  dB, 

PADB  =  Acoustic  power  per  primary  frequency  +0.5  dB. 
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The  PA  which  is  the  acoustic  power  per  primary  frequency  In  watts  Is  also 
proportionately  more  accurate,  as  the  +0.5  dB  tolerance  for  PADB  implies. 

The  new  formula  for  NDIS  as  well  as  the  formulas  for  BWS  and  BWP  depend 
on  GL  and  x  where 

GL  =  antilog  (G/20) 

and  x  is  selected  from  a  table  of  x  values  formed  by  the  program  according 
to  an  equation.  Because  GL  and  x  always  appear  as  the  ratio  x/(6L)»  the 
tolerances  on  NDIS,  BWS  and  BWP  will  never  exceed  +1.0  dB.  The  tolerances 
on  NDIS,  BWS  and  BWP  will  vary  and  depend  on  the  operating  point  on  the 
Mellen  and  Moffett*  curves. 

An  equation  for  determlng  a  table  of  x  values  has  been  mentioned.  This 
equation  determines  the  Increment  of  x  in  a  table  formed  by  the  program  and 

stored  in  the  computer  for  use  in  determining  the  values  of  LSS  through  the 

use  of  a  numerical  Integration  algorithm  by  Goldstein.4  The  table  Is  formed 
such  that  each  next  increasing  value  of  x  is  a  fixed  percentage  of  the 

previous  value.  In  order  to  reduce  the  tolerance  on  LSS  it  was  necessary 

to  reduce  the  fixed  percentage  Increase  on  successive  values  of  x*  The 
required  new  formula  for  x  is 
• 

X(I+1)  =  1,04881  xi 

for  172  values  starting  from  x0  B  0.00953462. 

New  Plots 


When  It  comes  to  analyzing  the  results  for  a  given  set  of  design 
parameters,  the  single  on-line  plots  of  LSS  vs  F0/F  and  NDIS  vs  F0/F  for 
each  different  combination  of  secondary  frequency  (F),  area  (A)  and 
nominal  source  level  (LSSI)  are  cumbersome.  What  has  been  devised  to  over¬ 
come  this  shortcoming  Is  plots  of  LSS  vs  F0/F  and  NDIS  vs  F0/F  for  a  given 
family  of  LSS  values  for  each  combination  of  F  and  A.  Figures  1  through  6 
show  the  results  of  this  new  plotting  method  for  the  sample  example  of 
Appendix  A.  The  plots  weie  made  on  an  off-line  plotter  (Calcomp)  that  is 
used  with  the  Uni vac  1108  computer. 

These  plots  can  save  many  man  hours  of  hand  plotting  to  achieve  a 
similar  matrix.  Also,  the  family  of  curves  approach  places  all  of  the 
data  for  a  given  F  and  A  on  one  piece  of  paper  and  allows  the  designer  to 
conveniently  view  all  the  available  trade-offs.  How  the  new  plotting  method 
fits  Into  the  computer  program  flow  Is  shown  In  Figure  7. 
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CONCLUSION 

A  computer  program  to  aid  In  the  design  of  parametric  sonars  has  been 
upgraded.  The  latest  Information  has  been  Incorporated  In  calculating  the 
secondary  frequency  directivity  Index  (NOIS).  Also  two  new  quantities  have 
been  Introduced.  They  are  the  secondary  frequency  3  dB  beamwldth  (BWS)  and 
the  primary  frequency  3  dB  beamwldth  (BWP).  These  three  quantities  are  valid 
for  a  circular  and  square  piston  projector. 

The  tolerances  on  the  calculated  Information  have  been  tightened.  Thus, 
for  example,  the  secondary  source  level  (LSS)  can  now  be  computed  to  +0.5  dB. 

A  new  plotting  method  using  the  family  of  curves  approach  Is  now  used. 
Plots  of  acoustic  power  per  tone  vs  downshift  ratio  (PADB  vs  FO/F)  and 
secondary  directivity  Index  vs  downshift  ratio  (NDIS  vs  FO/F)  can  now  be 
made  for  a  family  of  LSS  values  for  a  fixed  secondary  frequency  (F)  and  a 
projector  area  (A).  This  new  type  of  plot  enables  the  designer  to  more 
readily  observe  trends  and  to  make  tradeoffs. 
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FIGURES  1  THROUGH  3 
DESIGN  CURVES  (PADB  VS  FO/F) 
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FIGURES  4  THROUGH  6 
DESIGN  CURVES  (NDIS  VS  FO/F) 
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This  appendix  shows  a  sample  example  and  the  form  the  tables  and 
on-line  plots  take  as  they  come  out  of  the  computer.  The  off-line  plots 
were  shown  In  Figures  1  through  6  In  the  main  text. 
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